Andreev current in finite sized carbon nanotubes 
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We investigate the effect of interactions on Andreev current at a normal-superconducting junction 
when the normal phase is a Luttinger liquid with repulsive interactions. In particular, we study the 
system of a finite sized carbon nanotube placed between one metallic and one superconducting lead. 
We show that interactions and finite size effects give rise to significant deviations from the standard 
picture of Andreev current at a normal- superconductor junction in the nearly perfect Andreev limit. 



PACS numbers; 7I.10.Pm, 74.80.Fp, 73.63.Fg 

In recent years, the behavior of superconductors (SC) 
in contact with Luttinger liquids (LLs) has commanded 
attention in both theory [Q|2j and experiment [^,1). 
Josephson junctions made by sandwiching a Luttinger 
liquid between two superconductors have led to intrigu- 
ing results such as critical currents orders of magnitude 
larger than expected |^ . Experimental study of Andreev 
physics at a niobium superconductor - carbon nanotube 
junction [Q has yielded significant deviation at low tem- 
peratures from the standard picture of Andreev current 
in a non-interacting one-dimensional electron gas - super- 
conductor junction As it has been predicted Q and 
shown that single-walled metallic carbon nanotubes 
(NT) exhibit Luttinger liquid behavior, systematic anal- 
yses of any set-up involving their electronic properties 
would require taking into account the effect of interac- 
tions. 

Here we study the Andreev physics in a SC-NT- 
metallic lead junction, focusing on the effects of the 
strong repulsive interactions and of the finite size of the 
nanotube. We focus on energy scales well below the en- 
ergy gap A of the superconductor. Thus, throughout the 
energy range of interest (i.e for all the values of temper- 
ature or of the applied voltages we consider), the only 
excitations allowed to exit or enter the superconductor 
are Cooper pairs and not single electrons. In particular 
we focus on the limit of almost perfect Andreev reflection 
at the SC-NT interface (i.e. very low normal backscatter- 
ing). We also assume that the nanotube-metal contact is 
ideal and that the nanotube continues adiabatically into 
the metallic lead. Under these assumptions we study how 
a small amount of backscattering at the SC-NT interface 
would influence the electrical properties of the junction, 
in particular the behavior of the conductance. 

The treatment we use to obtain the value of the current 
as a function of the applied voltage is a non-equilibrium 
Keldysh technique, perturbative in the bare backscatter- 
ing strength 'm' [p|JTl|]. Characteristic of Luttinger liq- 
uids, the amount of backscattering can strongly increase 
when the energy at which the system is probed decreases. 
Hence, perturbation theory holds good only above an en- 



ergy scale Ec ~ eo(u/eo)^^^^^^\ where g measures the in- 
teraction strength { g = 1 in the absence of interactions). 
For metallic nanotubes, eo ~ \eV is the sub-band spacing 
and g ~ 0.25 |^J^, corresponding to strong repulsive 
interactions. In the set-up considered here, the effect of 
the finite length L of the nanotube becomes important 
below the finite size energy scale Tiv/L. Here, v = vp/g 
is the velocity of the charge-carrying quasiparticles in the 
nanotube, where vp ~ lO^m/s is the Fermi velocity. Ef- 
fects of finite size can be captured in the perturbative 
approach, as done here, provided Tiv/ L 3> Ec- 

To summarize our results, a numerical analysis re- 
veals that at zero temperature, the conductance shows 
a marked drop with decreasing voltage as a consequence 
of LL physics, consistent with renormalization group ar- 
guments similar to the ones derived in At voltages 
much smaller than the finite size energy hv/L, the con- 
ductance levels off to a constant. In addition, it exhibits 
small spikes with a voltage spacing of t:Tiv/2L (about 
2 — 3meV for a nanotube of micron length) , reminiscent 
of resonance peaks from quasi-bound states for charge 
carriers confined within the length of the tube. 

We now present the explicit calculation yielding the 
conductance as a function of applied voltage for the SC- 
NT-metal system described above. The s-wave SC lies 
in the region a; < and we assume it to be ideally con- 
tacted to a finite size nanotube of length L in the region 
< X < L which continues adiabatically into a metallic 
lead for x > L. We model the system in the semi-infinite 
region x > 0, as a four channel LL with interaction pa- 
rameters appropriate for the nanotube up to a; = L, and 
appropriate for no interactions for a; > L. The bosonized 
Hamiltonian for this system in the absence of any normal 
backscattering is given by 
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For simplicity, we have set the constants h ^ e = = 1. 
Note that a = p±,cr± correspond to the four free sectors 
of the theory and are obtained by linear transformations 



1 



from the spin-channel indices (1|,1J,,21",2J,) Q. The re- 
lation between the bosonic fields Oia, 4'ia (« = 1/2, a =t 
/ I) and the original chiral right-/left-nioving electron 
fields "^iR/La is expressed through the bosonization pro- 
cedure via the transformation "^iR/La ~ e*^'^'"^^'"-'. In 
the nanotube region < a; < L where interactions are 
present in the net charge density p+, we take gp+{x) — 
g « .25, and gp^,cr± — 1- Also, in the non-interacting re- 
gion X > L, we take ga{x) = 1 for all a's. The velocities 
of the free modes are given by Va{x) — vp/gaix). The 
total charge density is ptot = "^dxOp^/ir. 

In the almost perfect Andreev limit, the electrons in- 
cident from the nanotube side on the SC-NT interface 
reflect back as holes with opposite spin: /j,(0) = 

*U/t(0) and 'f^RyliO) = ^Lyt^O)' ^^^ere i refers 
to the channel indices 1 and 2. In the bosonized lan- 
guage, these boundary conditions become <j)p± (0) = 
and (0) = 0. 

The weak normal backscattering at the SC-NT junc- 
tion can be modeled by modifying the Hamiltonian to 
H = Ho+H', with: 

^"i E E [^lRo.m^uo) + h.c.] 

i=l,2a=T,i 

= wcos[0p+(O)]cos[0,_(O)], (2) 

where the bosonized form takes into account the Andreev 
boundary conditions at the SC-NT interface. For simplic- 
ity, we choose not to include the backscattering processes 
where particles can flip their band index since these terms 
do not give rise to any new physics. 

Following Ref. [0, we integrate out the (p variables 
in the action, as well as the entire ^-dependence away 



from a; = 0. The resulting unperturbed imaginary time 
Euclidean action becomes 

^o-^i: Mi^p.(-")p + ^ E i-"ii^^p-K)^ (3) 

where [3 is inverse temperature. Here the imaginary 
time Fourier transforms for all fields 'A' are defined in 
the standard fashion A{LJn) — Jq dr' A{T')e^'^"'^ , w„ = 
27m//3. The spatial variations of the interaction param- 
eter g{x) and of the velocity v{x) are reflected by the 
fact that the effective interaction parameter g{uJn) is fre- 
quency dependent and has the form: 

5K) = ^4^,with,^(i-^)e-l-l. (4) 

Here T = 2L/v is the time it takes a charge-carrying 
quasiparticle with velocity v to bounce back and forth 
between the ends of the tube. The limits L,w„ ^ and 
L,uJn — > oo retrieve the expected form g{u!n) = 1 and 
g{wn) = 9 for a semi-infinite Fermi liquid and a semi- 
infinite nanotube respectively |p^ . 

Along the lines of Ref. jl^], we proceed to construct 
a real time Keldysh action. We introduce 9f fields run- 
ning over forward and backward paths in time. We define 
da = {0+ + 9;^)/2a,iid9a = O+-0-. The resuhing action 
is S — So + Si + S2 where Sq is the unperturbed action, 
Si describes the effect of the weak backscattering at the 
SC-NT interface, and S2 captures the effect of applying 
a chemical potential difference V = dta. Thus 



J 2tt 2 g[uj) g{-uj) 2T ' J n g{Lu) 

f doj , , m9 f du! ~ , , . 

+ J ^^coth— IVHI + J — 

Si=iuJ di{cos[0+ (t)] cos[d+_ (t)] - cos[d;^ (t)] cos[e;_ {t)]}, 
f 2uj 

S2= da; — [aiu;)ep+{-Lj) + Tj{Lj)ep+{~u;)], (5) 



where g{u!) is the analytically continued version of 
g{ujn) in Eq.(^) with |w„t| replaced by iuir. For all 
fields 'A\ we have used the real time Fourier transform 
convention A{uj) ~ J dtA(i)e*"*. The source field 77 al- 
lows for calculation of the current, I{t) = 29p+{t)/Ti = 
{i / 2ti) [5 S / 5'q{t)[r,={)) . Average quantities may be derived 
by taking expectation values with respect to the Keldysh 
generating functional, ^ = / V[e+\D[e-]e-^ . 

Using the above Keldysh action, and treating the 
backscattering to lowest non-vanishing order in pertur- 
bation, we find the expectation value of the current to 
be 



/ = 8— F 
h 



(6) 



where from here on, we reinsert factors of e, 7i and fcs. 

The first term is associated with the constant conduc- 
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tance Go = 8^ in the absence of backscattering. As 
expected, this ideal Andreev conductance of the finite 
size nanotube in the presence of a metallic lead is that of 
a four mode non- interacting ID electron gas ||l|,|o|. The 
backscattering current Is takes the form: 



27r \h 



dt sin 



2eVt 



JC(t)-C(0)] 



sin[i?(i)]. (7) 
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In the above equation, C{t) = J2a=p± ^a{t) and 
-^(0 = Ea=p±[^a(-i) - Ra{t)]/2. For each mode, 
Ca{t) and Ra(t) are the correlation and response func- 
tions respectively, with Ca{t) — {da{t)Oa{0))o and 
Ra{t) = —i{(^a{t)(^a{0))o- Thcir Fourier transforms are 
related by the fluctuation-dissipation theorem, CaiLo) — 
— coth (^ ^^'^rp )lm[Ra{uj)]. Here averages are with re- 
spect to the unperturbed action, and we have i?p+(cj) = 
—iTTg(uj)/uj and Rp-iuj) — —i-k/lo 

We now make a series expansion of gioJn) in Eq.(||), 
qjujn) ~ Qy^,„ctny"-. When analytically continued, this 



gives: 

oo 

5H=g^a„(— |)"e— . (8) 

n=0 ^ 

Substituting the above in Eq.(|^) and taking the deriva- 
tive with respect to the applied voltage gives the follow- 
ing reduction in the conductance due to backscattering 
Gb = 4^, in the limit T ^ 0: 



k=0 n=0 ""^ L I V u / ; J ^^-^ 



where the coefficients /3„ = 5(q;„/2)[(1 — g)/{\ + g)]", 
for n > 0, and /?o = g)/2. We have used the high en- 
ergy cut-off eo to evaluate the Fourier transforms Ca{t) 
and Ra{t). Here the terms involving nr correspond to 
physical processes of n bounces of the quasiparticles at 
the boundaries of the nanotube. Note that besides the 
weak backscattering at the SC-NT interface, the quasi- 
particles can also backscatter at the nanotube-metallic 
lead junction even in the absence of a barrier, solely as 
a result of the mismatch of the values of the net charge 
and velocity of the free modes in the nanotube and in the 
metal. 
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FIG. 1. The net differential conductance G = dl/dV of 
the nanotube set-up in units of /h as a function of applied 
voltage eV in units of h/r. The values of the parameters are 
g = .25, u = 300(ft/r) and eo = W^ih/r). a), b) and c) 
show the conductance for different ranges of applied voltage. 
The log-log plot of Gb in d) shows a constant at low voltages 
and a power law Gb ~ V^~^ on average for high voltages, as 
indicated by the dashed line. 



The most revealing analysis of Eq.(|^) comes from nu- 
merical evaluation. We restrict the infinite sums of Eq. ^ 
to a finite number of terms, and introduce an explicit high 
energy cutoff eo, in order to regulate singularities. We en- 
sure that errors coming from both truncations are negli- 
gible. In Fig.l, we plot the net differential conductance 
— Gq — Gb, as the one of experimental relevance. The 
conductance drops with decreasing voltage and levels off 
at voltages much smaller than H/t. To see why this might 
be expected, notice that at large voltages eV :s> h/r, the 
time scale at which the system is probed is much shorter 
than T, and the conductance roughly behaves as if the 
nanotube were semi-infinite. Characteristic of Luttinger 
liquid physics, it thus drops as Gb oc u'^{2geV/eo)^~^ on 
average, as shown in Fig. 1 d). This limiting behavior 
can be seen directly in Eq.(^ by taking t — > oo. At 
low voltages eV <g; Ti/t, the associated time scale is long 
enough to capture the effect of the metallic lead and of 
multiple backscattering events at its interface. With de- 
creasing voltage, the conductance ultimately levels off to 
a constant, as per Ohm's law for the metallic lead. This 
limit can be obtained in Eq.(^ by setting r = 0. 

A striking feature of the plots is the presence of spikes 
at probe values eV = nhn/T, with 'n' being an integer. 
As thcir magnitude is minuscule compared to the net 
variation in conductance, it would be difficult to mea- 
sure them in experiment. However, these resonances do 
exist, and are signatures of the quasi-bound states that 
one would expect within the nanotube region, given that 
here the interaction parameter and velocity are different 
from those of the metallic lead. 

As a variation of the above set-up, let us now replace 
the nanotube by a finite sized Luttinger liquid with only 
two transport channels (spin t / i)- Such a situation can 
be realized by using, for instance, an etched quantum 
wire. The corresponding free modes carry net charge 
p and spin a, and are linear combinations of the two 
spin species. Andreev boundary conditions at the su- 
perconducting junction require 0p(O) = O,0a-{O) = in 
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corresponding bosonized variables. Thus, the system can 
be effectively described by a single channel in the 6 vari- 
ables. This allows for us to study the particular situation 
where the velocity of the charge mode in the Luttinger 
liquid V — vpj g would equal the Fermi velocity v^p in the 
metallic lead, i.e v^p — vp/g- Hence, we can focus on the 
physics arising purely from the mismatch of the charge 
of the elementary excitations in the Luttinger liquid and 
the lead. This would not have been possible for the case 
of the nanotube as it is described by two modes moving 
at different velocities, and matching the velocity of one 
mode to the Fermi velocity of the metallic lead would 
cause a velocity mismatch in the other mode. We cal- 
culate the conductance as a function of applied voltage 
for this system in a manner completely analogous to the 
one described above for the nanotube. The major differ- 
ence here is that we have only one mode with effective 
interaction parameter 

^^""^ ^ l-(l-.)exp(-Kr|)- (10) 

The resulting conductance is plotted in Fig. 2. The 
magnitude of the resonances at eV = nhir/T spans a 
larger fraction of the net variation in conductance com- 
pared to the case of the nanotube. It is noteworthy that 
these resonances exist in spite of the fact that there is 
no mismatch of velocities of the free phonon modes in 
the Luttinger liquid and in the metal. As expected, the 
phonons rebound at the Luttinger liquid-metal interface 
solely due to the impedance mismatch in the charge sec- 
tor. 
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FIG. 2. The net conductance G of the two-mode Luttinger 
liquid set-up in units of e'^/Zi as a function of applied voltage 
eV in units of h/r. The values of the parameters are g — .25, 
u = lA{h/T) and eo = W^ih/r). 

To summarize, we have looked at how the standard 
picture for Andreev current through a superconductor- 
metallic wire junction gets altered in the Andreev limit 
in the presence of interaction and finite size effects. The 
Andreev conductance shows a reduction with decreasing 
voltage which finally levels off at the lowest voltages. Fi- 
nite size effects also give rise to resonances manifest as 
small spikes in the conductance. 



Finally, turning to experiment, while Luttinger liq- 
uid behavior in nanotubes contacted to normal leads 
has been analyzed in great detail 0, by no means has 
it been studied systematically in the presence of super- 
conducting leads. As seen here, one would certainly ex- 
pect Luttinger liquid effects to yield significant deviations 
from the standard picture of Andreev physics for non- 
interacting one-dimensional wires. Consistent with our 
assumptions, such experiments can be performed, for ex- 
ample, in supcrconductor-nanotube junctions, in which 
the superconducting gap energy is of the order of several 
meV, while for a nanotube of a few microns, the finite 
size energy is in the range of a meV. At temperatures of 
the order of lOOmK, thermal effects are expected to be 
negligible. These conditions are well within experimen- 
tal reach, and systematic analyses of such set-ups could 
potentially reveal rich physics arising from bringing Lut- 
tinger liquids in contact with superconductors. 
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